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We compute the amplitudes for non-leptonic annihilation decays of B mesons into two particles 
within the pQCD collinear approach. The end point divergences are regulated with the help of an 
infrared finite gluon propagator characterized by a non-perturbative dynamical gluon mass consistent 
with recent lattice simulations. The divergences at twist-3 are regulated by a dynamical quark mass. 
Our results fit quite well the existent data of B" — >■ D~ K'^ and _B" -^ D~*K'^ for the expected 
range of dynamical gluon masses. We also make predictions for the rare decays B'^ — >■ K^ K^ , 
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I. INTRODUCTION 



Perturbative QCD has been largely applied in the 
study of B meson decays, specially in non-leptonic de- 
cays which requires the calculation of hadronic matrix 
elements. In such processes the long distance interactions 
between the initial and final state mesons impose a dif- 
ficult obstacle in amplitudes calculations. In the pQCD 
approach the calculation of the matrix elements can be 
performed since the perturbative dynamics can be factor- 
ized from the non-perturbative dynamics by applying the 
hard exclusive scattering Brodsky and Lepage's approach 
[ij, |2|] . In this procedure the non-perturbative contribu- 
tion in exclusive hadronic processes are contained in the 
wave functions that describe the hadrons as bound states 
of quarks. The employment of pQCD is justified due to 
the large scale involved in the strong interactions, since 
the energy of gluons exchanged in such processes are of 
the order of the B meson mass Q . 

There are some problems in this approach. One of 
them is the occurrence of soft divergences for the van- 
ishing gluon momentum, as well as collinear divergences 
when the momentum of the gluon is parallel to the mo- 
mentum of a massless quarks (since the masses of quarks 
u,d,s are usually neglected). Double logarithmic di- 
vergences also occur with both soft and collinear diver- 
gences. Usually, these divergences can be absorbed into 
the meson wave functions and the factorization is pre- 
served. However in the amplitude calculation of B meson 
decays may appear endpoint divergences that can not ab- 
sorbed into the wave functions. This happens when there 
are non-factorisable contributions related to interactions 
with the spectator quark and for annihilation interac- 
tions, respectively at twist-3 and twist-2. These contri- 
butions, in general, are power suppressed and are negli- 
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gible in front of the factorisable contributions. However, 
it is important to have a reliable calculation of the anni- 
hilation diagrams since they can generate strong phases 
that are relevant in the determination of the CP viola- 
tion parameters, as well as in the evaluation of decays 
that occurs through pure annihilation processes. 

In order to deal with these divergences several ap- 
proaches have being proposed in recent years. In pQCD 
approach by H. N. Li et. al [^-iSi] a kx factorization 
is performed, where the transversal momentum of the 
quarks are taken into account. The large logarithms are 
then resumed resulting in Sudakov factors that suppress 
the soft dynamics. In the QCD factorization approach 
(QCDF) proposed by Beneke et. al |9l-[l3l. the perturba- 
tive QCD is used at some extent. The authors argue that 
the transition form factors, that enter in the calculation 
of B meson decays, cannot be obtained through pertur- 
bative methods because they are dominated by "soft" 
interactions, and the Sudakov factors are not enough 
to suppress them. On the other hand, it is considered 
that the non-factorisable contributions are dominated by 
hard gluon exchanges. Therefore, in QCDF the terms 
that are dependent on the transition form factors are 
parameterized as in the naive factorization, and the non- 
factorisable contributions are perturbatively calculated 
in the collinear approximation. The endpoint divergences 
in non-factorisable diagrams, are treated with a phe- 
nomenological and model dependent parameterization as 
in the following equation 



^ = ln^(l + pe'^), 0<p<l. (1) 

X Ah 



Although this computation is an ad hoc one it allow us 
to use this approach to estimate the contributions of the 
non-factorisable diagrams. 

One third possibility that has been successfully applied 
several times in order to deal with the divergences in the 
collinear p QCD ap proach is the use of an IR finite gluon 
propagator [13l4l8l | , which appears as a non-perturbative 
solution of the gluonic Schwinger- Dyson equations (SDE) 



[19J and has been confirmed by lattice simulation of pure 
QCD [201 ■ This gluon propagator is characterized by a 
dynamical gluon mass, which naturally provides a self- 
consistent and model independent calculation. This ap- 
proach works well at leading twist [16] , but at twist-3 
coUinear divergences may also appear. These divergences 
occur in the massless quark approximation. Of course 
such divergences may be cured introducing current quark 
masses, which lead to an unnatural strong contribution 
from the fermion propagators of the light u and d quarks 
at twist-3. Moreover for consistency with the gluon dress- 
ing we cannot neglect the fermion dressing that gives dy- 
namical masses of order Aqcd to the light quarks. 

In this work, we will study the regulation of the diver- 
gences occurring in the coUinear factorization of pQCD. 
We focus on pure annihilation B meson decays, in which 
the dependence on the regulation will be more evident. 
We calculate the branching ratios of several pure anni- 
hilation channels of i?^, B^ and i?^. Our main goal is 
to revisit and improve our previous calculations at twist- 
2 |16j . using an infrared gluon propagator in agreement 
with the most recent QCD lattice simulations and per- 
forming the twist-3 calculation. In order to deal with 
the coUinear divergences that appear at twist-3 we will 
use the dynamical masses of the light quarks. This ap- 
proach is mandatory if we want consistency with the 
use of dynamical gluon masses. We study the follow- 
ing decay channels into pseudo-scalar particles: B^ — >■ 
TT+TT", D±7r=F, B° -^ K+K-, DfK^ and B+ -^ DfK°; 
and also the decay into a pseudo-scalar and a vector par- 
ticle, B° -^ D*'f'K^, 5° -^ Dl^K^ and B+ -^ D*+K°. 
We will then compare our results with the available data. 
We have to mention that these decays are quite rare, and 
most of them are beyond the present experimental lim- 
its. It is expected however that they can be detected 
with the advent of LHCb (LHC, CERN), or even earlier 
by the experiment CDF (Tevatron, Fermilab). 

The distribution of the paper is the following: In Sec. 

II we present the basic expressions of the pQCD approach 
to calculate the amplitudes of B meson decays. In Sec. 

III we show the annihilation amplitudes for the different 
B decays. In Sec. IV we discuss few aspects of dynamical 
mass generation and the regulation of the divergences. 
Section V contains our numerical results and Section VI 
is devoted to our conclusions. 



where Gp is the Fermi constant, Vckm are Cabibbo- 
Kobayashi-Maskawa factors, Qi are the four-quarks oper- 
ators contributing to the decay and Ci(/z) are the Wilson 
coefficients . The amplitude of the decay i? — >■ / is then 
obtained from: 



A{B^f) = {f\H,s\B) 



^E^CKMa(M)(Q. 



(3) 



where (Qi) — {j\Qi\B) are the hadronic matrix elements 
between the initial and final states. 

To compute the amplitudes we use the pQCD formal- 
ism in the coUinear approximation, where the hadronic 
matrix elements are obtained through a convolution of 
the hard scattering kernel and the distribution ampli- 
tudes of the mesons involved in the process. In the case 
of a two-body non-leptonic decay B — )■ M1M2, the matrix 
element of the operator Qi is given by: 



{Q^ 



Ax dy dz T, (x, y, z) $Afi {x) $m2 iv) $b [z), (4) 



where x,y,z are the momentum fractions, ^m{x) 
are the light-cone distribution amplitudes for the 
quark-antiquark states of the mesons, which are non- 
perturbative functions of the momentum fraction carried 
by the partons; Ti is the hard scattering kernel that can 
be perturbatively computed as function of the light-cone 
momenta of coUinear partons. The branching ratio is 
then calculated from 



B{B -)■ M1M2) 



tbPc 



\AiB -^ AfiM2)|^ 



(5) 



where tb is the B meson lifetime, Pc is the momentum 
of the final state particles with masses mi and 7712 in the 
B meson rest frame, 



1 



Pc = 



2m_B 



(m| 



(mi + rn2Y){m?B - (mi - 1112)'^)). 

_ , (6) 

As usual, it is used the light-cone distribution ampli- 
tudes decomposed in terms of spin structure. The de- 
composition of the wave functions of the B mesons are 
expressed in terms of two Lorentz scalar wave functions: 



II. PERTURBATIVE QCD APPROACH 

Weak decays of B mesons are described by an effec- 
tive Hamiltonian at a renormalization scale /x <?C Mw- 
This effective Hamiltonian is formally obtained using the 
Operator Product Expansion (OPE), describing effective 
interactions through four-quark operators, and is given 
by Hi: 






»(m), 



(2) 



'^B,{ali){x) 



--[{fl^)ali +"^-5(75)0/3] X 



X ( (t)B{x) + ^(I)b{x) ] , 



(7) 



where the normalization of the scalar functions are the 
following: 



dx(/)B(a;) =0. (8) 



dx(t)B{x) = — ^; 
2V6 



However, as it was shown in Refs. [22,|23|, the contribu- 
tion from the first wave function is dominant, therefore 




FIG. 1: Feynman diagrams contributing to the ampli- 
tude of anniiiilation non-leptonic decay channels. The non- 
factorisable diagrams are a) and b), and the factorisable ones 
are c) and d). 



a B meson is shown in Fig. [I] where the bottom quark 
decays through the W exchange or via a penguin process 
and a quark-antiquark pair is created in the final state 
through a ghion emission by any of the quarks involved 
in the process. The channels with these characteristics 
that we shall analyze are the following: 

• Charmless channels : 



• Charmed channels: 
B" ^ Dt^*'>K^; B^ 



^+(*)^o. 



In the following we present the analytical formulas for 
the computation of the amplitudes for the non-leptonic 
decays of B mesons. 



in our calculations we will consider only the first term of 
the Lorentz structure. 

For the D pseudoscalar mesons we use the following 
decomposition |23l |: 



^D{aP){x) 



[inbf)ali + "2d(75)q/j] (I)d{x) ■ (9) 



The light pseudo-scalars mesons are represented by the 
following wave function 



$ 



K,-iT[aP) 




- l)'/>T(a;)]a/3, 



(10) 



where the wave function (j)Aix) is the twist-2 wave func- 
tion and 4'p{x),(J)t{x) are twist-3 wave functions, and 
mop = Mp/{mqi -t- m,2) (toott = 1-5 GeV, mox = 
1.6 GeV). The vectors n,v are light-like vectors: for a 
meson moving in the direction P oc (1, 0, 0^), the vector 
n — {1,0, Ot) defines the direction of motion of the me- 
son, and we define v ~ (0, 1,0^). For a meson moving 
with momentum defined along the v direction, the third 
term of Eq. (|10[). we have that ^^ is replaced by ^^. 

For the heavy vector mesons (D't^^), only the longitu- 
dinal part of the wave function contributes to the ampli- 
tude, due to the conservation of angular momentum. We 
use the following expression: 



$? 



^pix) 



where e 



1 

V2r 



(1, 



= ^[(l^75)a^ +mD'{l5)af)](l}D'ix), 

(11) 

r^, Ot) is the longitudinal polarization 



of the vector meson, with r = mu/niB- 

III. ANNIHILATION AMPLITUDES 

Annihilation decays are processes that occur at as or- 
der in perturbative QCD. A typical annihilation decay of 



A. Charmless channels 

We present the amplitudes for the diagrams of Fig. [TJ 
obtained with the help of Eq. ^. The diagrams that 
actually contribute for the amplitude of charmless de- 
cays are the non-factorisable diagrams shown in Figs. [T^ 
and[TlD, since the contributions from the factorisable ones 
(Figs.[T]: and[T]i) cancel among themselves. The full am- 
plitude for the decays B —>■ tttt, KK is given by 



A{B -> M1M2) 



fsfMi fAh (Kb V* A 



^tb *^tq '^Penguin)-) 



Tree 



(12) 



with q = d and q = s for the decays of the -B^ and i3° 
mesons, respectively. In particular the amplitude for the 
decay mode i?^ — >■ tt^tt'^ relates to the one of the mode 

TT+TT- as A{B°s -^ 71'"'''") = 72"^(^° -^ Vr+TT-). 

We consider the B meson at the rest frame, with 
momentum Pb — ^^{1,1, Ot) in hght-cone coordi- 
nates (p+, p^,pt), with p"^ — -y={p° ± p^) and px — 

{p^,p^) (pt = Ot in the coUinear approximation). The 
mesons Afi and M2 have momenta defined as Pi — 
^(0,1, Ot) oc V and P2 = ^(1,0, Ot) oc n, hence 
Pb = Pi + P2- We also define the momenta of the 
quarks in the mesons B, Mi, M2 respectively as q = 
(zP+,0,Ot), gi = (0,xPf ,0t), 92 = (yP2+,0,OT) (the 
momentum of the antiquarks are given hy q — P — q). 
The tree level and penguin amplitudes are 



A' 



•^Pcne 



2C4 



Tree 

C'lO 
~2~ 



= C2A 



Al 



2Ce 



2 



(13) 



A2. (14) 



The functions Ai and A2 are the contributions from 
operators of the type {V - A) (» {V - A) and {V - A) (» 
(V + A), respectively, and are given by 



Ai 



=" Ann- / dxdydz(pB{z)<hh(x,y,z} X 

\JlNc Jo I 

x(j)A{x)(l)A{y) + Mop ( -'^<Pp{x)(j)p{y) + {z - y) x 



D,[*', K, n) will be calculated from 



AiB -^ D-M^ 



Gp 
Gi 



Vu,V:^M{B^D-M+), 



A{B^D+M-) ^ -^VcgV^k M{B ^ D+M~) {19) 
X {(t>p{x)+4)T{x){(t>p{y)-(t)T{y))+x{(t>p(x)~4)T{x))x v2 



X {(f'piv) + 4>T{y)) 



+ hq{x,y,z) 



+ f^lp S(0p(x) + cl>Tix)){My) - ^T{y)) + 



+ (y - z){M^) - 0T(a;))(<^p(2/) + Mv)) 

and 

32TTm%GF f\ , , , , J , , 

A2 ^ _— — / dxdydz(pB(zj<.hb(x,y,z)x 

^/zlMc Jo L 

X i-y + z-l) (t)A{x)(j)A{y) + IJ-lp { -A(t)p{x)(j)p{y) + 

+ x{4)p{x) + c^T{x)){4>p{y) - My)) + {z-y) {<t>p{x) + 



{y — z)(j)A{x)(j)A{y) + where q — diox the decays of B^ and i3+, and q — s for 

We consider the B meson at the rest frame, with mo- 
mentum Pb = ^(1, 1, Ot). The light mesons M ^tt,K 
1 /ir\ have momenta defined as Pm — ^(0, 1 — ?'^,0t) oc v, 

} and for the charmed mesons the momentum is Pd = 



Pd 

2i|(l, r^, Ot), hence Pr = Pm + Pd- The momentum of 
the quarks are defined analogously to the charmless case. 
The functions M in the previous Eqs. ([T^ are the 
sum of factorisable and non-factorisable contributions, 
M = Miac+ -A^nfac- The analytical formulas for these 
contributions for the decays with two pseudo-scalar par- 
ticles in the final state, are the following: 



0T(x))(0p(y)+0T(y)) 



+ hq{x,y,z) 



X(j)A{x)(t)A{y) 



• B" 



-M+- 



+Mop f (y - z){(i))p{x) + ^T{x)){(t>p{y) - My)) + 

+x{(j)p{x) - (pT{x)){4>p{y) + 4>T{y)) 

Where fiQp = mop/mB, and the functions h{x,y,z) are 

hb{x,y,z) = a{fj.)Di,{x,y,z)Dg{x,y) 
hq{x,y,z) = a{nh)Dq{x,y,z)Dg{x,y), (17) 

and the quarks and gluon perturbative propagators are 
defined as: 



D-M 



> ^.-.,^+^_ 8TTm%CF '^ 



(16) Xfac(i?° ^ D-M+) = 



dx dy dz 



V2iVc Jo 
X (j)B{z)(l)D{y){ -( {1 - r'^)y(l>A{x) ~ 2rfiop{l - r'^ + y)(l)p{x) 



X hi{x, y, z)- f (1 - r'^){r'^ - (1 - r'^)x)(j)A{x) - r^Lopi^ - r"^ + 
- 2(l-r2)x)0p(x)-(l-r2)(l-2x)0T(a^)Vc(x,y,z)l, (: 



Db^{x,y,z) = m%x{z-y)-ml 
Dq^{x,y,z) = m%x{y-z)-ml 



D-\x,y) = (fc^)-i 



m%yx 



(18) 



The scales adopted in the calculations are /i — nib for the 
contribution from the diagram shown in Fig. [T^ where 
the gluon is emitted from the h quark, and ^u = V^Uijnb 
for the contribution shown in Fig. [TJj where the gluon is 
emitted from the spectator quark, with Ah = 500 MeV 
Q. When the effective scale is fi < mb we use Aqcd — 
225 MeV and when /i = nib we use Aqcd = 300 MeV, 
and the difference is to match the values of the coupling 
constant due to different quark thresholds. 



B. Charmed channels 

The charmed decay channels occur through W ex- 
change processes, and there is no contribution from pen- 
guin diagrams. The full amplitude for the general pro- 
cesses B -^ D^M^ (with B = Bd,B°,B+; D,M = 



MniUB°^D-M+)^^^^^^^ f dxdydzM^)My: 



c Jo 



X |- ((1 - '^')((1 + r^){y -z)- l)Mx) - r/iop((2 + y + z- 

+ (1 - r^)x)(l)p{x) + {-y + z + {l- r^)x)(l)T{x)) j hb{x, y, z) 4 

+ ( (1 - r'^)iiy - z)r'^ + (1 - r'^)x)(t)A{x) + rfiQp{{y -z + 

+ (1 - r^)x)4>p{x) + {y-z-{l- r^)x)(f>T{x)) ] hg{x, y, z) \ .( 



• B° ^ D+M- 



MiUB" -^ D+M- 



-MfaciB'^ -^ D-M+) (22) 



Xnfac(S" ^ D* + M- 



M.,..{B^^D+M-)=''^<''^ '' 



32TTmgCF 

V2]Vc JO 



dx dy dz x 



, , dx dy dz x 

V2iVc Jo 

X ,^B(z)0c(y)|- ("(1 - r2)((l - r^)x + {y- z^^Aix) + 

- r^0p((2 + 2/ + z + (1 - r2)x)0p(.T) + ((1 - r^)x + 

+ y + z)4>t{x)) ) hbix, y, z)+ ( (1 - r^){y - z)4>a{x) + 

+ r^J.op{{y - z + (1 - r^)x)(f)p{x) + 

+ {-y + z + {l- r^)x)(l>T{x))jhg{x, y, z) \ 

For decays with one vector plus one pseudo-scalar par- 
ticles in the final state, the analytical formulas are the 
following: 

m B° ^ D*-M+: 



(23) 



x(t>B{z) My) 1 - ( (1 - ^')(i -iy + ^y - (1 - '^')-^) ^ 

X(f)Aix) +rnop{{y + z - (1 - r'^)x)(l)p{x) - {y + z + 
-2r2 - (1 - r2)(l + x))Mx))\hi{x. y, z) + 

+ m - 2r2)(y - z)x)(I)a{x) + r^op((2/ - z - {I - r'^)x) x 

X(/)p(a;) - (y - z + (1 - r'')x)<i>T{x))\hq{x, y, z)\ (27) 



In the charged decay modes -B+ — )• Ds -^° the an- 
alytical formulas for the amplitudes are easily obtained 
from the ones of the neutral modes i? -^ Dg K~, 
given in Eqs. (20), (21), (24) and (25). The simpler way 
to see this, is that the operators that contributes to the 
neutral decay modes are related to the ones that con- 
tributes to the charged modes by replacing the spectator 
quark of the B meson d ^- u (i.e., bd — > 5m), and the 
light quark from the K meson is also replaced u ^ d 
(i.e., us — >■ ds). The final result in the calculation of 
„i the amplitudes of the charged modes are equivalent to 

/ dxdydz(f)B{z) x the neutral mode, with the replacement of the distribu- 
"'0 tion, <j>BO ~^ 4'B+ Slid (t>K -^ 4'K°- Also the analytical 

, / X J / /, 2\ ± f \ n /I 2 \ J. / \\ expression of the distribution amplitudes are equivalent, 

X 4>D[y)< {{l~r )y(j)A[x) ~2rfioP[^-r - yy)(j)p[x) ]x 



MtUB"" -^ D*-M+) = 



%/2]Vc 



X hi{x,y, zy (1 - r2)(l - (1 - r'')x)(f>Aix) 



+ r^op((l - r'))<l)p{x) - (1 + r')(t)Tix) Kix, y, z) 



with the use of the respective mass and decay constant 
of each meson. Besides that, the amplitudes of decay 
modes which involves two neutral pseudoscalar mesons 
in the final states are related to the modes of charged 
j'24'i pseudoscalars as 

A{B -^ M^M°) = -^A{B -> M+M^) 
v2 



■Mnfac(i3° ^ D*-M+) = ^^"^gl^"^ /' dxdydz X 

V2iVc Jo 

X <l>B{z)My)l-(i^ - '■')((! - r^)iy + z) + r^)cbA{^) + 

-rpiop{{y + z-{l- r^)x)(f>p{x) + {y + z ~ 2r^ + 

- {l-r^){l+x))(l)Tix))jhb{x,y,z) + 

+ ({l-r^)iil-r^)~r^{y~z)x)Mx) + 

- rnop{{y- z- {l~r^)x)(j)p{x) + 

+ (y- z + {^-r^)x)(jJTix))]hg{x,y,z) 



In Eqs. (PH)) to ([?7)) we define the following functions: 

hc{x,y,z) = Dc{x) Dg{x,y)as{iJ,f){Ci{fif) + C2{fif)/^) ; 
hi{x, y, z) = Di{x)Dg{x, y)as{fif){Ci{^if) + C2{fif)/3) ; 
hg{x,y,z) = Dg{x)Dg{x,y)as{^ih)C2{^J.h) ; 
hb{x, y, z) ^ Db{x)Dg{x, y)as{p)C2{n) , (28) 

where the functions Di,{x,y, z), Dd{x,y, z), Du{y) and 
Dc = Dc{x) and Dg{x, y) are the quarks and gluon prop- 
agators, given by: 



Dy{x,y,z) = ml{{y - -z){l - x{l - r^))) 



'6 7 



(25) Dq i^iy^z) = mpx{z-y){l-r ); 



'))-m. 



2 . 
c 5 



• B" 



D*+M- 



D-^(x) = m|(r2-hx(l 
DiHy) = mly{\-r^) 
D-\x,y) - kl = mly{l-x){l-r^) 



(29) 



For the charged modes, the Wilson coefficients that ap- 
pear in Eqs. (28) are exchanged by Ci(///)/3 + C2(/i/) 
A^fac(B° -^ D*+Ar) = -Mfac{B° -^ D*-M+) (26) for the first two hues, and by Cl{^l) for the last two hues. 



Note that the perturbative propagators are going to 
be substituted by the non-perturbative ones discussed 
in the next section. The scales adopted in the calcula- 
tions are the same as the charmless case: fi — mb for the 
contribution of the diagram where the gluon is emitted 
from the b quark, and fih = \/ ^h'nih for the contribution 
where the gluon is emitted from the "spectator" quark, 
with A/i = 500 MeV. For the diagrams of Fig. [U; and[TJi 
where the gluon is emitted from the quarks in the final 
state, we use the scale /x/ = /z/2 = mb/2. The QCD 
scales are also the same as the charmless case. 



IV. REGULATION OF THE DIVERGENCES 

The annihilation amplitudes have both soft and 
collinear divergences. The soft divergences occur with 
the vanishing of the gluon momentum, and this one has 
been treated with an ad hoc cutoff as shown in Eq.([T|). In 
our work this procedure is exchanged by the introduction 
of an IR finite gluon propagator obtained through the 
solution of Schwinger-Dyson equations (SDE) [19l . y2A - 
u7\ and consistent with recent lattice data [20] • As for 
the collinear divergences, we can include the dynamical 
mass of the light quarks, which are also necessary to keep 
the consistency with the introduction of dynamical gluon 
masses. The poles of the (massive) propagators will give 
the imaginary part of the amplitudes. Therefore, the 
basic point in the regulation of divergences is based on 
the use of non-perturbative IR information (or dressed 
propagators) in the context of the perturbative expan- 
sion. We will discuss this regulation scheme in detail in 
the following. 



A. IR finite gluon propagator 

A prescription of how non-perturbative SDE solutions 
can be inserted into the perturbative QCD expansion 
was proposed by Pagels and Stokar many years ago, in 
the approach denominated dynamical perturbation the- 
ory (DPT) [2g]. In their scheme the amplitudes that do 
not vanish to all orders in perturbation theory are given 
by their free field values, while amplitudes that vanish 

1/2 

as A cx e~ '^ are retained, and possibly dealt with in 
an expansion in g"A. The work of Ref.[28] was partic- 
ularly concerned with the effect of a dynamically gener- 
ated quark mass, but from SDE [1^, [24l - [27f and lattice 
simulations [20] we now know that the gluon and cou- 
pling constant also have an infrared finite value, and the 
Pagels and Stokar formulation can be extended and gen- 
eralized for the perturbative expansion using the quark 
and gluon propagators with a dynamical mass and the 
IR finite coupling constant. Actually we expect that any 
infrared finite gluon propagator leads to a freezing of the 
infrared coupling constant [29], meaning that the use of 
an infrared finite gluon propagator must be accompanied 
by an IR finite coupling constant. All these facts indi- 



cate that we should perform perturbation theory with 
the dressed quark and gluon propagators and the effec- 
tive charge (dependent on the gluon mass). 

The SDE solution and the lattice result that we shall 
use were discussed at length in the following references 
|20l [30| , therefore we will not enter into further details 
about the solutions and will just present the gluon prop- 
agator and coupling constant in the case that QCD gen- 
erates a dynamical gluon mass. We consider a gluon 
propagator that will have the form 



zA^,(g)=P^,A(<z)+^M^ 



Pn 



-9,.+^ , (30) 



where A(g) is the gauge invariant scalar part of the gluon 
propagator, which in Euclidean space has the form 



A(Q") (X 



1 



r+m^Q^) 



(31) 



The gluonic SDE solutions allow us to write a new prop- 
agator A~^((5^) which absorbs all the renormalization 
group logs, exactly as happens in QED with the pho- 
ton self-energy, and form the product d{Q^) — g^A{Q^) 
which is a renormalization group invariant. The dynam- 
ical gluon mass {m^{Q^)) [1^ is given by 



ml{Q^) = ml 



In 



(^^) 



In 



4m2 



(32) 



where the gluon mass scale nig has typical values [13, |311 



mg = 500 ± 200 MeV 



(33) 



As this is a complicated expression to take into account 
when performing the numerical integration, in practical 
calculations, the best we can do is to assume rni{Q^) sa 
rri^. In this work the dynamical gluon mass is calculated 
at the same scales described in the end of the Section 
III. There is no significant difference in the results if we 
neglect the running of the masses, since slightly away 
from the mass value the effect of momentum dependence 
takes over the full propagator effect" . 

A simple fit for the coupling constant that is factored 
out in this procedure is given by [191 ] 



asdiq"^ 



1 



47r61n[(4m2 



q"^ - ie)/K^ 



(34) 



where 6 = (33 - 2n/)/487r2. Eq.dMD clearly shows the 
existence of the IR fixed-point as g^ — > 0. It must be 
stressed that the fixed point does not depend on a spe- 
cific process, it is uniquely obtained as we fix A and, in 
principle, it should be exactly determined if we knew how 
to solve QCD. Several examples of the use of DPT with 
the propagator and coupling constant discussed above 
can be found in Ref . Isil . 



B. Dynamical quark masses 

The dressed quark propagator can be written as 



S{p) = 



t-f ■p + M{p'^) 



(35) 



where Z{p^) is the Lorentz scalar quark wave function 
renormahzation and M{p'^) is the quark dynamical mass. 
At large momenta, and in the chiral limit, it is known 
that that the dynamical mass is related to the quark 
condensate {< qq >) as [S^l 



Af(p2 



<qq> 

p2 



(36) 



The many estimates of the dynamical quark mass (Af (0)) 
(among them we have M(0) (x < qq >'^^^) give values be- 
tween 250 MeV and 300 MeV ^ [11], whereas we can 
safely assume Z{p^) « 1. The decrease of the mass with 
the momentum shown in Ea. (l36p is not going to be con- 
sidered, since the mass function has a plateau at low 
momenta and decreases (in the chiral limit) very fast 
at high momenta 33], where the propagator is domi- 
nated by the ij ■ p term, in such a way that our cal- 
culation is not affected if we take the quark propagator 
as S{p) = l/{ij-p + M{0)). 



The Wilson's coefficients are computed using the equa- 
tions given in the appendices of Ref. [g . We also use the 
following parameters |36| : 

• Masses: 

TUB = 5.28 GeV, niB, = 5.37 GeV, mo = 1.87 
GeV, niD, = 1.97 GeV, niD- = 2.01 GeV, mo* = 
2.10 GeV, TO6 = 4.7 GeV. m^ = 1.27 GeV. 

• Lifetimes: 

tb^ = 1.54 ps, tb^ = 1.466 ps, tb+ = 1.638 ps . 

• GKM parameters: 

A = 0.814, A == 0.2257, p = 0.135, 77 = 0.349. 

• Decay constants [3J, |33 : 

/s = 190 MeV, /b, = 236 MeV ; fo = 206 MeV, 
fa, = 257 MeV, /^ = 130 MeV, Jk = 156 MeV; 
fa- = 245 MeV, /r,; = 272 MeV. 

In Table U we show the results obtained and the exper- 
imental data available for each decay channel. 



VI. CONCLUSION 



V. NUMERICAL RESULTS 

We perform the numerical calculation considering the 
gluon mass scale with values m.g — 400 MeV, 500 MeV. 
The dynamical masses of the light quarks are M (0) — 
niu = rud, in the interval 250 MeV < Af (0) < 300 MeV, 
and the same value is assumed for the s quark, although 
it has been argued that the constituent s quark mass may 
be slightly higher |33|. For the the light mesons, pions 
and kaons, we use the asymptotic expression for the wave 
functions 



4F{x) 



t P P t P 

6x (1 — cc), (p {x) 



2V2iVc 
fp 



2V2Nc 



{2x - 1) 



2V2^ 
(37) 



while for the D{D*) mesons we use 



'^d{x) = 



h 



2y/2N, 



=6.t(1-.t) [1 + 0^(1 -2a;)], (38) 



c 



with a£i — 0.8 for D and a£i — 0.3 for Dg mesons 
|38l I40J . The momentum of the B meson is carried al- 
most all by the b quark, and the light quark carries a 
momentum of the order Aqcd- Providing that we as- 
sume z ^ ^Qco/'mb and x,y ^ z, at leading order the 
integration over the B meson wave function is trivial, 
yielding to the decay constant according to the Eq. ([5|). 
This procedure has been adopted in annihilation decay 
calculations by Beneke et. al |9|]. 



We have studied some decay channels of 5°, 5°, B^ 
mesons which occur through the annihilation diagrams. 
We have argued that infrared finite gluon propagators 
and running coupling constants obtained as solutions of 
the QGD Schwinger-Dyson equations, may serve as a nat- 
ural cutoff for the end-point divergences that appear in 
the calculation of these decays. Also, we use the con- 
stituent quark mass for the light quarks in the propaga- 
tors to be able to perform the twist-3 calculation of the 
amplitudes. 

Comparing to the data available, our calculations are 
very satisfactory, notably for the two channels with ac- 
tual measurements B'^ — >• Dj K+ and B° — >• Dj*K+. 
The results are also reasonably compatible with some 
predictions of the pCJCD in the kr factorization, see for 
example Refs. |37h43]. 

We can see from the results that the variation on the 
dynamical mass of the quarks generally introduces only a 
small error. The numerical results are more dependent on 
the gluon mass scale, as we can see a more significant vari- 
ation of the results with this mass scale. Comparing to 
the available data it is possible to narrow the interval for 
the gluon mass scale for values closer to rUg ~ 500 MeV. 
This value for the gluon mass scale is in accordance 
with previous calculations [31j . Future measurements 
for the other decay modes, besides i?" — > D^K'^ and 
B'^ — > D~*K~^, can establish a more precise value for the 
nig parameter. Besides the annihilation modes, it is also 
possible to apply the present regulation prescriptions to 
other types of B decay. 



TABLE I: Branching ratios for B decays. The error correspond 
•niq < 300 MeV. We show the results for the gluon mass scale nig 
available experimental data. 



to the variation of the dynamical quark mass, 250 MeV < 
= 400 MeV and 500 MeV. In the last column are shown the 



Channel i 


ng = 400 MeV 




mg = 500 MeV 




Experimental 
















Data ;30] 


B" ^ K+K- 


;4.31±0.50) 


X 


io-« 


(3.20 ±0.01) 


X 


io-» 


< 4.1 X 10^' 


m -^ 7r+^- 


;i0.42±1.00) ; 


<io-^ 


(5.32 ±0.25) 


X 


lO"'^ 


< 1.7 X lO-** 


Bl -> 7r%° 


;5.21±0.05) 


X 


10-^ 


(2.66 ±0.13) 


X 


10-^ 


< 2.1 X 10-'' 


B" ^ D-K+ 


;6.12±0.78) 


X 


10"^ 


(3.07 ±0.17) 


X 


10-5 


2.9 ±0.5 


B° -> DtK- 


;6.15±0.71) 


X 


io-« 


(2.92 ±0.23) 


X 


io-» 




B+ -> D+RO 


;0.89±0.16) 


X 


io-« 


(0.32 ±0.11) 


X 


io-» 


< 9 X 10-'' 


fiO ^ D-7T + 


;i.72±0.04) 


X 


lO"*^ 


(1.18 ±0.40) 


X 


lO"'' 




B° ^ 5%o 


;0.86 ± 0.02) 


X 


lO"*^ 


(0.61 ±0.17) 


X 


lo-'' 




fiO ^ D+7T- 


;5.80 ± 0.60) 


X 


10"^ 


(2.51 ±0.96) 


X 


lO"'^ 




B° ^ D%o 


;2.90±0.30) 


X 


10-^ 


(1.26 ±0.48) 


X 


10-^ 




B° ^ Dl~K+ 


;2.85±0.18) 


X 


10"^ 


(1.90 ±0.06) 


X 


10^5 


2.2 ±0.6 


B° ^ Dl+K- 


;6.23 ± 0.50) 


X 


10"* 


(2.68 ±0.07) 


X 


10"* 




B+ -> Dl+KO 


;i.71±0.52) 


X 


io-« 


(1.28 ±0.22) 


X 


io-» 


< 9 X 10"^ 


B'l -^ D*-Tx+ 


;3.36 ± 0.06) 


X 


lO"*^ 


(1.88 ±0.21) 


X 


lO"'' 




Bl ^ D*+7r- 


;3.26±0.29) 


X 


io-« 


(3.92 ±1.31) 


X 


io-» 





The already existent agreement of some B decays with 
the experimental data and the possible future confirma- 
tion of the other rare decays computed here, may in- 
dicate that the introduction of the nonperturbative in- 
formation provided by the Schwinger- Dyson equations in 
the form of dynamical masses for quarks and gluons, con- 
sistent with the most recent lattice QCD simulations, is 
gathering more and more robust phenomenological re- 
sults showing that these mass scales possibly cannot be 



neglected when performing high precision hadronic cal- 
culations. 
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